Energy decay laws in strongly anisotropic MHD turbulence 
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We investigate the influence of a uniform magnetic field Bo — -Boe|| on energy decay laws in 
incompressible magnetohydrodynamic (MHD) turbulence. The nonlinear transfer reduction along 
Bo is included in a model that distinguishes parallel and perpendicular directions, following a 
phenomenology a la Kraichnan. We predict a slowing down of the energy decay due to anisotropy in 
the limit of strong Bo, with distinct power laws for energy decay of shear- and pseudo-AIfven waves. 
Numerical results from the kinetic equations of Alfven wave turbulence recover these predictions, 
and MHD numerical results clearly tend to follow them in the lowest perpendicular planes. 
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Predicting the evolution of freely decaying turbulence 
remains one of the most difficult problems in turbulence. 
For the sake of simplicity, this problem is often tack- 
led under the assumption of homogeneity and isotropy. 
In the case of 3D Navicr-Stokcs flows, the Kolmogorov 
prediction for the kinetic energy decay law thus reads 
Evit) ^ {t — t*)~^°/'', where is the time origin for 
power law decay Close power law indices arc mea- 
sured in grid turbulence experiments although boundary 
effects may alter the result 0]. The generalization of 
decay laws to other situations is still currently under de- 
bate. For example in rotating turbulence, a recent exper- 
iment shows that neutral fluids behave differently from 
pure (non rotating) flows with a slowing down of the ki- 
netic energy decay due to the presence of anisotropy 3 . 
In this case, the original isotropic Kolmogorov descrip- 
tion has to be modified to include rotation effects that 
may lead to strong spectral anisotropy (see e.g. Ref. [3]). 

In this Letter, we investigate the influence of an exter- 
nal uniform magnetic field on the energy decay laws in 
freely incompressible magnetohydrodynamic (MHD) tur- 
bulence. The MHD approximation has proved to be quite 
successful in the study of a variety of astrophysical plas- 
mas like those found in the solar corona, the interplan- 
etary medium or in the interstellar clouds. These me- 
dia are characterized by extremely large Reynolds num- 
bers (up to 1013) d with a range of available scales from 
lO^^m to few meters. The isotropy assumption is partic- 
ularly difficult to justify when dealing with astrophysical 
flows since a large scale magnetic field is almost always 
present like in the inner interplanetary medium where the 
magnetic field lines form an Archimedean spiral near the 
equatorial plane (see e.g. Ref. @). The present study, 
although theoretical, appears therefore particularly im- 
portant to extract some universal features of turbulent 
plasmas. The MHD equations in presence of an external 
uniform magnetic field Bo= -BoG|| read 

atv-Bo9|tb + v- Vv= -VP, + b-Vb + t^Av, (1) 
9tb-Bo9||V + v-Vb = b- Vv + ryAb, (2) 



with V • V = and V • b = 0. The magnetic field b is nor- 
malized to a velocity (b — > y/^Qnrrii b, with the ion 
mass and n the electron density) , v is the plasma flow ve- 
locity, the total (magnetic plus kinetic) pressure, v the 
viscosity and r] the magnetic diffusivity. The role of the 
Bg-field on the flow behavior has been widely discussed in 
the community (see e.g. Ref. 0, i, B 0, El, El, El, H ) ■ 
One of the most clearly established results is the bi- 
dimcnsionalization of an initial isotropic energy spectrum 
with a strong reduction of nonlinear transfers along Bq. 

In the past, several papers have been devoted to predic- 
tions of energy decay laws in isotropic MHD turbulence 
[H, El, El, Ei, El- We first review one of them [l| 
which follows the Kolmogorov-Kraichnan spirit that we 
then adapt to strongly anisotropic flows. In the simplest 
case (balance turbulence), the derivation of self-similar 
decay laws in MHD relies on few basic assumptions which 
are (i) a weak correlation between velocity and magnetic 
fields (which allows to use below the variable z instead 
of , where = v ± b), (ii) a power law spectrum 
in E{k) ~ k'^ (with s = D -\- 1, D being the space di- 
mension) for the low wavenumber energy, i.e. at scales 
larger than the integral scale i from which the inertial 
range begins, and, obviously (iii) a power law time de- 
pendence for E{t) ~ {t — t^)~"' and l{t) ~ (t — t^,)^, where 
a and /3 arc two unknown indices. Another hypothesis, 
directly related to the invariance of the Loitsianskii inte- 
gral [201, t'^Us that the modal spectrum scales like A:""^ 
at low wavenumbers and, with sufficient scale separation, 
this dominates the total energy so that E ^ ^~(^+i\ 
hence the first relation a = (3{s + 1). A second re- 
lation may be obtained from the energy transfer equa- 
tion e = —dE/dt ^ E/Tfr ~ (i — t*)~"~-^, where Ttr is 
the transfer time and e is the transfer rate. According 
to the Iroshnikov-Kraichnan (IK) phenomenology [2l| . 
Ttr ~ ^j^hl'^^^ with TATi = Ij zi, the eddy turnover time, 
and TA ~ ^/Po, the Alfven time. Substituting previous 
times into the energy transfer equation leads to 1 = q;+/3. 
Finally one obtains a = (s + l)/(s + 2) and (i — l/(s + 2). 
The predictions for 3D MHD turbulence (s = 4) are then 
E(t) ~ (t-t,)-5/6 and^(t) - (i-t^fl^. The energy de- 
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cay law for isotropic MHD turbulence has been favorably 
compared with direct numerical simulations, in particu- 
lar in the 2D case for which higher Reynolds numbers 
can be reached than in the 3D case [l^, [l^ . Note that 
other phenomenologies exist for explaining, in particular, 
the influence of magnetic helicity on the decay laws [l^] , 
a situation different from here where it is negligible. 

We now derive an anisotropic version of the previ- 
ous heuristic description which is well adapted to MHD 
flows permeated by a strong uniform magnetic field. The 
main characteristic of such flows is that nonlinear trans- 
fers are strongly damped along the mean field Bq direc- 
tion, leading preferentially to perpendicular (_L) trans- 
fers. In Fourier space, this means that the modal en- 
ergy spectrum mainly develops a power law scaling in 
A;^-wavenumbers whereas the scaling along parallel wave 
numbers, fcy, does not change very much. In practice, 
we assume that the intensity Bq is strong enough to ig- 
nore parallel transfers and we make the approximation 
k ^ k± fc|| which means, in particular, that the initial 
condition of the system is confined to the largest parallel 
scales. In this situation, it is straightforward to show that 
shear-Alfven wave energy corresponds to perpendicular 
fields (z_l) whereas pseudo-Alfvcn wave energy mainly 
comes from parallel fields (z||). Our heuristic description 
is focused on shear-Alfven waves, and not pseudo-Alfvcn 
waves, since it is well known that under such anisotropic 
configuration, the latter are slaved to the former [l^l . 

We make three equivalent hypotheses as in the 
isotropic model and we note E±, £± and £|| the shear- 
Alfven wave energy, the perpendicular and parallel in- 
tegral scales respectively. The length scale, although 
appearing sometimes in our derivation, is assumed time 
independent because of the strong nonlinear anisotropy. 
In that case, we assume (a) v± ~ 6^ ~ (b) for fc|| 
fixed, E±{k±, fcy) ~ k^ (where D is the space dimension) 
at low perpendicular wavenumbers, i.e. at scale larger 
than e±, and (c) E±{t) ~ (t - and £_i^ {t- t*)^, 

where a and (3 are the new unknown indices. The sec- 
ond assumption means a, D — 1 power law index for the 
modal spectrum before the integral fcj^-wavenumber; fur- 



thermore it leads to estimate E±{t) 
hence the (same) first relation 

a = P{D + l). 



[D+l] 



in , and 



(3) 



A second relation may be obtained by using the energy 
transfer equation 



= -dE^/dt - E^/Ttr - (t - t*)" 



(4) 



where Ttr = 't'^l/'''A: with, for anisotropic transfers, 
tnl = ^/ik±z±£^) and ta ~ l/(fc||i?o)- Substituting 
these times into Q leads to the (new) relation 

1 = a + 2/3 . (5) 

Finally, ^ and ^ lead to two new scaling exponents: 



£±1 -R 1 



Hence the predictions for 3D anisotropic MHD turbu- 
lence {D = 3) E^[t) - (t-t,)-2/3 and£j_(t) - {t-U^/^. 
These results show, in particular, a slowing down of the 
energy decay for shear-Alfven waves compared to the to- 
tal energy in the isotropic case (where no distinction is 
made between shear- and pseudo-Alfven waves). 

The extension of our approach to pseudo-Alfven waves 
is not direct since these waves are (mainly) slaved to 
shear-Alfven waves. We remind that shear-Alfven and 
pseudo-Alfven waves are the two kinds of linear pertur- 
bations about the equilibrium, the latter being the in- 
compressible limit of slow magnetosonic waves. For the 
pseudo-waves, the heuristic description based on nonlin- 
ear transfers is misleading. Instead, it seems suitable to 
find a relationship between the different type of waves. 
The divergence free condition provides this relation which 
eventually leads to a prediction for the energy decay law. 
In Fourier space, the divergence free condition reads 







(7) 



where and Z|| are the Fourier transform of the 
cartesian fields which may be associated, under strong 
anisotropy assumption, mainly to the shear- and pseudo- 
Alfvcn waves respectively. Since we assume a weak cross- 
correlation (balance turbulence), it is not necessary to 
introduce the Elsasser variables. Simple manipulations 
lead to k'^E± ~ ^jj^Hi where iSy denotes the pseudo- 
Alfven wave energy. Since nonlinear transfers along the 
Bg-direction are negligible, may be seen as a mute 
variable. Therefore, the energy decay law for pseudo- 
Alfven waves should be 
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(8) 



13 + 3 ' D + 3 
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Our anisotropic model thus predicts quite different en- 
ergy decay laws for shear- and pseudo-Alfven waves, the 
latter not depending on the system dimensions. An im- 
portant issue concerns the balance ^ fc^^"^ law found in 
many simulations (see e.g. psl). Our analysis is based on 
the permanence of big eddies [l| which is directly linked 
to the conservation of the spectral scaling law at the 
largest scales, i.e. at scales larger than the ones of the in- 
crtial range where the "critical balance" is well observed. 
The 2/3 law is a priori not included in the model (as well 
as the exact scaling law for the energy spectrum) because 
the assumption of the existence of an inertial range where 
energy is evacuated from the reservoir is enough. There- 
fore, the integral length scales £± and are not linked by 
the 2 /3 law even in the derivation of equation ([5]) since we 
are only dealing with scales before the inertial range. The 
possible persistence of the 2/3 law at the largest scales is 
nevertheless an important issue which, in principle, may 
be investigated from the derivation of the von Karman- 
Howarth equation in anisotropic MHD turbulence, from 
which a Loitsianskii type invariant could be found. This 
is a huge challenge for the future which is, of course, out 
of the scoop of this paper. 

In order to check the model validity, we perform nu- 
merical simulations of strongly anisotropic MHD flows. 
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The first data set comes from integrations of the kinetic 
equations of Alfven wave turbulence derived and ana- 
lyzed in ■ This regime describes the asymptotic limit 
of strong Bq. We do not re-derive or even re- write these 
equations that the reader can find in [l^ (equations (54) 
and (55)). One of the main results found in this regime 
is a kj^ scaling law for the energy spectrum (for both 
shear- and pseudo- Alfven waves) which is an exact power 
law solution of the wave kinetic equations (in absence of 
cross-correlation). Another important result is the to- 
tal absence of parallel transfer. For that reason, the 
numerical simulations are made at a fixed fcju A non- 
uniform grid in Fourier space (see for details [IJl) is used 
to achieve a highly turbulent state at a Reynolds number 
of about 10^, with v = -q = 2x 10~^. Figure[T]displays the 
time evolution of the shear- Alfven, pseudo-Alfvcn and 
total energies in absence of cross-correlation. We first 
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FIG. 2: Temporal evolution of shear-Alfven energy spectra. 
A k\ and power laws are ploted for comparison. 




FIG. 1: Temporal decays of shear-Alfven (long-dash), pseudo- 
Alfven (dash-dot) and total (solid) energies, together with 
three theoretical slope predictions: 3D isotropic (3Di) in 
^-5/6^ 2D isotropic (2Di) in t""''^ and 3D anisotropic (3Da) 
in t^'^^^ . A t^^ slope is also ploted for comparison with the 
pseudo- Alfven energy decay. 

note the existence of a transient period during which the 
energy is conserved. During this period the energy cas- 

— 7/3 

cades towards smaller scales following a fc_|_ scaling law 
rather than the exact solution [l^. A clear power 
law behavior appears for both types of waves. These 
power laws are in agreement with the theoretical predic- 
tions made for 3D anisotropic flows (for comparison, the 
predicted slopes for the 3D and 2D isotropic cases are 
also given). Surprisingly, the total energy follows very 
precisely the theoretical prediction over more than two 
decades whereas pseudo- Alfven waves decay slower than 
at very large times. This discrepancy may be linked 
to the saturation of the integral length scale as it can be 
observed from the temporal evolution of the shear-Alfven 
wave spectra (Fig. [2]). Indeed, a self-similar energy de- 
cay is observed with a slow increase of the integral length 
scale that can be roughly estimated from the wavenum- 



ber at which the inertial range begins (i.e. the maxima 
of the spectra). At later times, this scale is close to the 
maximum size of the numerical box whereas the large- 
scale power law in k\ is still preserved. The large scales 
are even more reduced for (slaved) pseudo- Alfven waves 
since they decay faster. Actually, the finite size box ef- 
fect may explain the change of decay law at very large 
times for pseudo-Alfven wave energy. Direct numerical 




F+(k,, = 0) 



run 512^x64 Bn=15 



FIG. 3: Evolution of shear- (solid lines) and pseudo-Alfven 
(dash lines) wave spectral energies at fixed fcy wavenumber, 
together with the total energies. Note that curves at fcy = 
and fc|| = 1 are shifted by a factor of 2 for clarity. (Color 
online) . 



simulations of 3D incompressible MHD equations UHl]) 
are also performed with a pseudo-spectral code includ- 
ing de-aliasing. A high resolution with 512^ in the Bo- 
transverse planes whereas only 64 grid points are taken in 
the longitudinal direction. Such a situation was analyzed 
to explore the self-consistency of the reduced MHD model 
[23 | with the conclusion that small values of viscosities, 
adjusted according to the transverse dynamics, are not 



4 



incompatible with the smaller spatial resolution in the 
longitudinal direction since the transfer towards small 
scales is also reduced along the uniform magnetic field. 
We have checked that viscosity values, v = rj = 5 x 10^"*, 
are indeed well adjusted 25[. The initial condition cor- 
responds to a modal energy spectrum in agreement with 
the phenomenology described above, with a D = 3 power 
law at largest scales : E^{k±,k^^) = C{k^^)k^ for fcj^ and 
€ [0,4], the value of C(fc||) increasing with fcj| to reach 
a maximum at fcy = 4. This initial spectrum allows a 
transient period of cascade towards smaller scales during 
which energy is mainly conserved. The uniform magnetic 
field is fixed to Bq = 15. Initially, the ratio between ki- 
netic and magnetic energies is one, whereas the correla- 
tion 2(v • b)/(|vp + |b^|) is zero (remaining less than 4% 
up to t = 40). In Fig. O time evolutions of shear- and 
pseudo-Alfven wave energies are ploted at fixed parallel 
wavenumbers, namely fcy = 0, 1, 2, 3, together with to- 
tal energies (as obtained by integrating over all 
Note that we chose to show the behavior of one type of 
waves, i.e. y (with by definition = (|z^p)/2), as 

II energies behave the same way. We observe that sys- 
tematically E^ > E'^ at the final time which confirms 
the previous results. A clear tendency towards power 
law behavior is also found, in particular at low k\\ , with 
exponents close to the heuristic predictions. This was 
already observed at smaller resolution (not shown) and 
seems to be a general tendency of anisotropic MHD flows. 
However, this is less clear in planes at fcii =2 and 3, prob- 



ably due to enhanced dissipations there (for example, at 
t = 20, 60% of the shear wave energy is lost at fcy = 2, 
and ^ 73% at fc|| = 3) leading thus to shortened self- 
similar decay ranges. Actually, this remark may explain 
the power law steepening of the energy decay, when in- 
tegrated over all k\\ , since the energy loss is even more 
pronounced in higher fc|| -planes. This average effect has 
never been emphasized in the literature but seems to be 
the most important obstacle to see the decay laws as well 
as the spectral la ws p redicted by anisotropic theories like 
wave turbulence [2y|. Note that other initial conditions 
with an extended energy spectrum in fcy might reduce the 
discrepancy found between the model and the simulation 
at high /c|| , a situation not investigated here because of 
the constraint due to the reduced parallel resolution. 

In this Letter, we show the influence of a uniform mag- 
netic field on energetic decays in MHD turbulence. Mod- 
ified self-similar laws are derived, with a t~''^^^ and a 
decay, respectively for shear- and pseudo-Alfven waves, 
under strong anisotropy. To our knowledge, this decay 
analysis is the first in the context of wave turbulence and 
could be extended to many other problems. Integrations 
of kinetic equations of Alfven wave turbulence recover 
the predicted laws, while our MHD numerical flows fol- 
low them in planes at lowest parallel wavenumbers. 
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